Shielding of a small charged particle in weakly 



ionized plasmas 

Manis Chaudhuri, Sergey A. Khrapak, Roman Kompaneets, and Gregor E. Mornll 



Abstract — In this paper we present a concise overview of 
our recent results concerning the electric potential distribution 
around a small charged particle in weakly ionized plasmas. A 
number of different effects which influence plasma screening 
properties are considered. Some consequences of the results are 
discussed, mostly in the context of complex (dusty) plasmas. 

Index Terms — Complex (dusty) plasmas, screening, potential 
distribution. 



I. Introduction 

Complex (dusty) plasmas are plasmas containing small 
charged particles of solid matter (dust grains). These particles 
are usually large enough to be observed individually, which al- 
lows experimental investigation with high temporal and spatial 
resolution. Hence, complex plasma is a valuable model system 
for studying various phenomena (e.g, phase transitions, self- 
organizations, waves, transport, etc.) at the most elementary 
kinetic level 0, 0, 0, £Q, 0. 

The character of the plasma-particle as well as the inter- 
particle interactions appears to be one of the most fundamental 
questions for understanding the physics behind the observed 
phenomena in laboratory complex plasmas 0, 0, 0, 0, 
as well as in astrophysical plasmas [6|, plasma of fusion 
devices Q, 0, @, plasma processing 0, (TO), flj], etc. 
Of particular importance are basic processes such as parti- 
cle charging, electric potential distribution around a charged 
particle in plasmas, interparticle interactions, momentum and 
energy transfer between different complex plasma components 
etc. Ifl2l . 

The main aim of this paper is to present a brief overview 
of our recent results (mostly theoretical) related to the electric 
potential distribution around a charged test particle in plasmas. 
A number of different effects which influence plasma screen- 
ing properties are considered, such as plasma absorption on 
the particle surface, non-linearity of ion-particle interaction, 
ion-neutral collisions, plasma production and loss processes 
in the vicinity of the particle. These effects influence and 
often completely determine the shape of the electric potential 
around the particle, especially its long-range asymptote. This 
is important for a number of collective properties of complex 
plasmas (e.g., interparticle coupling, phase transitions, phase 
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diagrams, transport, etc). Below, we first discuss the role of 
these effects in isotropic plasma conditions (Section ITXb and 



then consider the case of anisotropic plasmas (Section III » 



This is followed by a short summary in Section IV 



II. Isotropic plasmas 

The distribution of the electric potential around a small 
individual spherical non-absorbing particle of radius a and 
charge Q in isotropic plasmas is often described by the Debye- 
Hiickel (Yukawa) form 

cj>{r) ~ <f) s (a/r) exp [-(r - a)/\ D ] ~ (Q/r) exp(-r/A D ), 

(1) 

where Ad is the linearized Debye length, Ad = 
Adi/vI + (^Di/^De) 2 - The ion (electron) Debye length is 
defined as ADi( e ) = \/Ti( e ) / Aire 2 ^^, where Tjf e ) is the ion 
(electron) temperature, ~ n e ~ n is the unperturbed plasma 
density, and e is the elementary charge. The particle surface 
potential <j) s = (Q/a) exp(— a/Ao) is usually of the order of 
the electron temperature, (f> s ~ —T e /e, due to much higher 
electron mobility. The above expression can be obtained by 
solving the linearized Poisson equation with the assumption 
that ion and electron densities follow Boltzmann distributions 
and the condition |e0 s /2^( e ) | < 1 is satisfied. It is to be 
noted that usually in complex plasmas T e 3> Tj and the 
linearization is often invalid since ion-particle coupling is very 
strong close to the particle. Nevertheless, numerical solution 
of the non-linear Poisson-Boltzmann equation shows that the 
functional form of Eq. ([TJ still persists, but the actual value of 
the particle charge should be replaced by an effective charge 
which is somewhat smaller than the actual one in the absolute 
magnitude fljl, lfl4t 

The effect of plasma absorption on the particle surface 
strongly modifies the potential distribution. The continuous 
ion and electron fluxes from the bulk plasma to the particle 
make their distribution functions anisotropic in the velocity 
space. Although the deviations are negligible for repelled 
electrons 031 . for attracted ions they are quite substantial. 
In the absence of plasma production and loss in the vicinity 
of the particle, conservation of the plasma flux completely 
determines the long-range asymptote of the potential and in 
collisionless plasmas it scales as <f>(r) cx r~ 2 0, lfl5ll . |16|. 
Close to the particle (up to a distance of few Debye radii from 
its surface), the Debye-Hiickel (DH) form works reasonably 
well. However, in the regime of strong ion-grain coupling, 
the linearized Debye length Ad should be replaced by the 
effective screening length A c ff. The exact dependence of A e ff 
on plasma and grain parameters is not known for the general 
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case and so far it was determined using numerical simulations 
only for a limited number of special cases [17|, [18], [i9J. For 
a collisionless plasma with T e ^> Tj a fit based on numerical 
results of Ref. ifTTl has been recently proposed [ 12]. The 
corresponding expression is A e ff — Ad[1+0.105\//5+0.013/3], 
where f3 = |Q|e/(TjAD) is the so-called scattering parameter 
[20 1, which is a natural measure of nonlinearity in ion-grain 
interaction. This fit has been used to calculate the ion drag 
force in a collisionless plasma with strong ion-grain coupling 
in Ref. ED and to estimate the effect of polarization inter- 
action on the propagation of dust acoustic waves in complex 
plasmas in Ref. [22] . 

Another important factor which influences the structure 
of the electric potential around an absorbing particle in 
plasmas is ion-neutral collisions. In the weakly collisional 
regime (£i > Ad, where £j is the ion mean free path) the 
infrequent ion-neutral collisions create trapped ions which 
increase the ion density in the vicinity of the particle. This 
affects the ion flux collected by the particle El], ED, E21 
and, therefore, modifies the potential distribution. In strongly 
collisional plasmas (£i <C Ad) the electric potential is known 
to exhibit a Coulomb-like (oc r" 1 ) decay |25), ED, ED, ED- 
A transition from the DH to the unscreened Coulomb potential 
with increasing ion collisionality was observed in a numerical 
simulation l30l . 

Recently, a simple linear kinetic model has been proposed 
independently by Filippov et al. Oil and Khrapak et al. Q"2l 
which takes into account the combined effect of ion absorption 
on the particle and ion-neutral collisions. Using this model 
the electric potential distribution can be calculated in the 
entire range of ion collisionality. Below we present the general 
expression for the potential obtained in Ref. 0"2l and analyze 
some interesting limiting cases. 

A. Linear kinetic model 

In this model a small (point-like) negatively charged indi- 
vidual grain immersed in a stationary isotropic weakly ionized 
plasma is considered. Plasma production and loss processes 
are neglected in the vicinity of the particle except on the 
particle surface, which is fully absorbing. The electron density 
satisfies the Boltzmann relation, and ions are described by 
the kinetic equation accounting for ion-neutral collisions and 
ion loss on the particle surface. Collisions are modeled by 
the Bhtanagar-Gross-Krook (BGK) collision integral with a 
velocity independent (constant) effective ion-neutral collision 
frequency v. Ion loss is expressed through an effective (veloc- 
ity dependent) collection cross section a. The linear response 
formalism is used to solve the Poisson equation along with the 
corresponding equations for ions and electrons. The resulting 
expression for the electric potential is 



(r) = — cxp(-fc D r)-- / 
r r J 



fc D sin(kr)f(9)dk 
fc 2 + fc 2 



9l + <PlI, 

(2) 



where 



/(*) 



8n / °° q(QC 2 arctan(C/fl) exp(-C 2 R 
7r 3 / 2 fc D 1 - ^9 exp((9 2 )[l - erf (9)} ' 



Here fc D i( e ) 
length, fco 
thermal velocity. 



— A Di(e) 



= vw 



is the inverse ion (electron) Debye 
is the ion 



C 2 = v 2 /2v 2 



rrii is the ion mass, 9 = (u/y^kv^), and 
■j t . . The first term fa in Eq. d2| is the familiar 
Debye-Hiickel potential. The second term <pn appears due to 
ion absorption by the particle and accounts for ion-neutral 
collisions. For a non-absorbing particle [a = 0] only the 
conventional DH form survives, as expected. In this case ion- 
neutral collisions do not affect the potential distribution. 
Let us now consider some important limiting cases. 
1) Collisionless limit: In the collisionless (CL) limit 
one can use the OML collection cross section a(Q = 
ira 2 [1 + [zt)C 2 ] to get EH 



<M r ) : 

where J-(x) = [e~ x 'Ei(x) 



to get 

e 7ra 2 n(l + 2zt) 



2k T 



T(kvr), 



(3) 



— e x Ei(— x)] and Ei(x) is the expo- 
nential integral. Here, z = \Q\e/aT e is the normalized particle 
charge and r = T e /Ti is the electron-to-ion temperature ratio. 
For sufficiently large distances x 1, T{x) « 2/x and the 
corresponding potential is 



Mr) 



T 



1 + 

4(1 



2zt 



(4) 



which coincides with the well known result of probe the- 
ory 0, 031, HSl, OS. The long-range asymptote of the 
potential scales as oc r~ 2 in this collisionless limit. 

2) Strongly collisional limit: In the opposite strongly col- 
lisional (SC) regime the actual form of is not important 
since the integral in f(9) is directly expressed through the ion 
flux Ji in this case. The resulting potential is 

e Ji 

<M0 ^ — [1 - exp(-fc D r)] , (5) 

where Di = v^./u is the diffusion coefficient of the ions. 
This expression coincides with the results obtained using the 
hydrodynamic approach (34], 031 . The long-range asymptote 
of the potential decays as cx r^ 1 in this case. 

3) Weakly collisional limit: The most interesting regime 
relevant to many complex plasma experiments in gas dis- 
charges is the weakly collisional (WC) regime, li > Ad- In 
this case the functional form of a(Q is required to calculate 
the potential distribution. A simple assumption of a constant 
cross section cr(£) = cr = \Z"7r/8( Jj/nwyJ was made in 
Ref. [32], which allowed to avoid divergence of the integrals 
in calculating f(9). The result is 



Ji 



r 4^2 knVTi 



^ + Sd- [1 



exp(-fc D r)] 



(6) 



The two terms in the curly brackets of Eq. (|6]l correspond to 
absorption induced "collisionless" and "collisional" contribu- 
tions, respectively. The collisional contribution to the potential 
dominates over the collisionless one for r > 0.6£i. 

Figure [T] demonstrates that the long-range asymptote of the 
potential is dominated by the combined effect of collisions 
and absorption. It exhibits Coulomb-like decay (f>{r) ~ Q e s/r, 
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Fig. 1. Distribution of the normalized electric potential around a small 
individual charged particle in an isotropic weakly ionized plasma for different 
values of the ion collisionality index. The solid curves are obtained using 
numerical integration of Eq. |5J. Dashed (dotted) curve corresponds to the 
analytic approximation in the strongly collisional (collisionless) limit. The 
dash-dotted curve shows the Debye-Hiickel potential with the surface potential 
calculated from the (collisionless) OML model. The inset shows a comparison 
between direct numerical integration of Eq. [2] (solid lines with symbols) and 
analytical approximation for the weakly collisional regime [Eqs. |2) and jSJ]. 



where the effective charge Q c g is determined by the plasma 
and particle parameters and increases monotonically in ab- 
solute magnitude with ion collisionality. At short distances 
the potential follows the DH form (QJ, but the actual particle 
charge Q shows a non-monotonic dependence on ^jfeo. In 
the WC regime \Q\ decreases with increasing collisionality, 
while in the SC regime \Q\ increases until it reaches a 
certain maximum value when both ion and electron collection 
by the particle becomes collision dominated ll36ll . Note that 
in the WC regime the transition from short-range DH to 
the long-range Coulomb-like asymptote can occur through 
an intermediate oc r -2 decay, whilst in the SC regime the 
potential is Coulomb-like practically from the particle surface. 

B. Effect of ionization/recombination 

So far it has been assumed that there are no plasma 
sources and sinks in the vicinity of the particle except at its 
surface. Physically, this corresponds to the situation when the 
characteristic ionization/recombination length is considerably 
larger than the characteristic size of the plasma perturbation 
by the charged particle, i.e., compensation of plasma losses to 
the particle occurs very far from it. However, in real conditions 
some plasma production and loss processes always operate in 
the vicinity of the particle and therefore, it would be interesting 
to estimate the importance of these effects. 

For an individual particle this has been done in Refs. ||29l , 
[37 1 using the hydrodynamic approach for the case of highly 



collisional plasmas. The effects of plasma production and loss 
processes are introduced in the ion continuity equation through 
the corresponding plasma source and loss terms. Electron 
impact ionization is usually considered as the main mechanism 
of plasma production. Plasma loss can be either due to 
electron-ion volume recombination ll29l (which is relevant to 
high pressure plasmas) or due to ambipolar diffusion towards 
the discharge chamber walls and electrodes [37] (which occurs 
in low and moderate-pressure gas discharges). In addition, 
when many particles are present in the system, plasma losses 
can be associated with absorption on the particles themselves. 
This situation has been recently analyzed in detail in Ref. [38|. 
Let us consider these three cases separately. 

1 ) Losses due to volume recombination: The plasma source 
and loss terms added to the ion continuity equation are in 
this case Vjn e — v^n e rii, where Vj and vr are the effective 
ionization and recombination constants, respectively. In the 
unperturbed state vi = v^n. The standard linearization pro- 
cedure yields the electric potential of the form 



(j> = (Q+/r) cxp(-rk + ) + (Q_/r) cxp(-rfc_), (7) 



where 



kl = 




and 



Q± ==F- 



Q [k% - k D 2 - (eJi/QDi)] 



k\ 



k 2 



(8) 



(9) 



The potential is screened exponentially but unlike in Debye- 
Hiickel theory it is described by the superposition of the two 
exponentials with different inverse screening lengths k + and 
k-. Both these screening lengths depend on the strength of 
plasma production (ionization frequency Vi). The effective 
charges Q + and Q_ also depend on plasma production 
strength as well as on the ion flux collected by the grain. 
The long range asymptote of the potential is determined by 
the smaller screening constant fc_ with effective charge Q_. 
Depending on the strength of plasma production two limiting 
cases can be considered: low and high ionization rate. 

In the limit of low ionization rate, vi/Di <C fc^, the screen- 
ing is dominated by ionization/recombination effects and the 
screening length Ado(^^d)\/ 1 V 1 ^ i s considerably larger than 
the electron Debye length since (likjf) 3> \Jvjjv in the con- 
sidered regime. For distances Ad < r < Xr> e (£ikr>) yvjvi 
the potential behaves as Coulomb-like with the effective 
charge Q_ ~ —{eJij Dik\f), i.e., we recover the result of the 
previous (no ionization/recombination) limit [Eq. (|5j]. Thus, 
the distance \v e (£ikv) yjvjvj determines the length scale 
below which plasma production is not important and sets up 
the upper limit of applicability of the results obtained within 
the assumption of no ionization/recombination processes in the 
vicinity of the grain [Eq. Q]. 

In the opposite limit of high ionization rate, vi/Di 3> fc^,, 
the screening length is given by the electron Debye length 
Aoe and is independent of the ionization rate V[ . The effective 
charge Q_ ~ Q — eJi/vj is somewhat larger in the absolute 
magnitude than the actual charge. 
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2) Losses due to ambipolar diffusion: In this case the 
plasma source and loss terms are vjn e — vi,ni, where vl is the 
characteristic loss frequency {vj — v£). The potential around 
the grain is l37l . 



4> = (Qi/r) exp(-fc eff r) + Q 2 /r, 



(10) 



where Q 1 = Q[l - {v T - (e/Q)J i }/D l k 2 s ], Q 2 = Q-Qi = 
(Qvi — eJi) I Dik^ff, and k 2 s = k^ 2 + vj/Di. Thus, in the 
considered case the potential is not completely screened, but 
has Coulomb-like long-range asymptote. The effective charge 
Q2 depends both on the strength of ionization vi and the ion 
flux Ji collected by the grain. The analysis of the limits of 
low and high ionization rates is straightforward. 

3) Losses due to absorption on other grains: In this case, 
the role of ion sink is played by a continuous, immovable, and 
uniform "particle medium". The ion source and loss terms are 
essentially the same as in the previous case, with a minor 
difference that the ion loss frequency vl may depend on 
the particle charge which varies with the ion and electron 
densities. The principal difference from the previous case is 
that n,i 7^ n e in the unperturbed state due to the presence of the 
charged particle medium. This has been shown to change the 
Coulomb term Q 2 /r in Eq. (IOi to (Q%/r) cos(kor) as well 

l39l. RUl The parameter k 



as to change Qi 2 and k e s (T9 
is , 



ko 



Pzv L 



rDi[P + z + {dv L /dz){z/v L )\ 



(11) 



where P = (rij — n e )/rii is the Havnes parameter (see, e.g.. 



Eq. (27) of Ref. ROT ). Equation (11 1 is derived for v L jD l < 



> 1, 



1, (dvL,/dz)(z/vi,) ~ 1 and yields ko 



10 2 &d for typical experimental parameters (see also Fig. 3 
of Ref. Ull). 

However, it has been recently shown that this cosine-like 
potential cannot be observed in principle (38). The reason 
is that such an ionization-absorption balanced plasma with 
rii ^ n e is unstable with respect to ion perturbations and 
the threshold wavenumber is exactly equal to ko, whereas 
the attraction was derived by implicitly assuming the ion 
component to be in a stable equilibrium and considering the 
screening of a test particle as a static perturbation of this state. 

This instability disappears for a constant ionization source 
[38 1 . However, for a constant ionization source the cosine in 
the potential changes to the exponent |41|, |42|. 

III. Anisotropic plasmas 

Electric fields are often present in plasmas (e.g., in rf 
sheaths, positive column of dc discharge, dc discharge stri- 
ations, ambipolar electric field in plasma bulk, etc.). This 
induces an ion drift and, hence, creates a perturbed region of 
plasma density downstream from the particle - the so-called 
"plasma wake". One can apply the linear dielectric response 
formalism 1431 to calculate the potential distribution in the 
wake. This approach is applicable provided ions are weakly 
coupled to the particle (i.e. the region of nonlinear ion-particle 
electric interaction is small compared to the plasma screening 
length). Note that higher ion drift velocities imply better appli- 
cability of the linear theory. The electrostatic potential created 



by a point-like charge at rest is defined in this approximation 

as 



(r) 



dk 



Q_ 

2tt 2 J fc 2 e(0,k)' 



(12) 



where e(uj, k) is the plasma dielectric function. Using a 
certain model for the dielectric function, one can calcu- 
late the anisotropic potential distribution using analytical ap- 
proach ED, 05J, 06), 03, 081, 09), (501. The potential 
profile can be also obtained from numerical simulations [51 1, 
1521 . Il53l, 1541 . 1531 . 11561 . 1571 . In general, the shape of the 
wake potential is sensitive to the ion flow velocity, plasma 
absorption on the particle, ion-neutral collisions, etc. Let us 
illustrate how the wake potential can depend on plasma con- 
ditions using few recent examples from analytic calculations. 

The examples below deal with a homogeneous plasma with 
ion flow driven by an electric field, where the unperturbed 
velocity distribution of ions is determined by the balance 
of the electric field and collisions with neutrals (mobility- 
limited flow). Here, "homogeneous plasma" presumes that the 
inhomogeneity length (e.g. due to Boltzmann distribution of 
electrons in the field driving the flow) is large enough. Such 
a plasma has been shown to be stable with respect to the 
formation of ion plasma waves, except for the parameter range 
where the thermal Mach number M? = u/vtu is larger than 
~ 8 and the ratio of the ion-neutral collision frequency to the 
ion plasma frequency, v /uj p i is less than ~ 0.2 [58|. Here u is 
the ion flow velocity, Vxn = \/T n /mi is the thermal velocity 
of neutrals, T n is the neutral temperature, and the ion plasma 
frequency is uj pi — y / 4ne 2 n/m i . This justifies the use of 
the linear response formalism ( 12 1 outside the aforementioned 
instability range. However, note that the stability analysis was 
performed using the BGK collision term, whereas the use of 
the more realistic constant-mean-free-path collision term 1591 
may yield somewhat different instability thresholds. 

A. Subthermal ion drifts 

In this subsection we focus on the subthermal ion drift 
regime, Mt < 1. 

The screening in the "almost collisionless" case has been 
investigated using the kinetic approach with the BGK collision 
term [60|. The "almost collisionless case" presumes that the 
ion-neutral collision length is much larger than the Debye 
length and distances where we want to find the potential. This 
case is treated by taking the limit v — > at a fixed Mt- In 
this limit the role of collisions is only to determine the non- 
Maxwellian form of the unperturbed velocity distributions of 
ions. The resulting potential for r S> Ad is l60l 



(r) = Q 



cxp(-r/A D ) /2M T A D 

— £\ I 7; — COS ( 



2 - 



(1 - 3 cos 2 



, (13) 



where 9 is the angle between r and the ion flow, and 
\e> = \jT n j '(Airne 2 ) (note a difference with respect to 
previous notation); the electron component is treated as a 
homogeneous neutralizing background which is not perturbed 
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Fig. 2. Distribution of the normalized electric potential behind a small 
individual charged particle in a weakly anisotropic plasma. The solid curve 
represents the potential distribution in the collisionless (CL) limit, Eq. j!3) , 
for \Q\e/(aTi) = 2 (here T, = T n ), M T = 0.2, and r > 1.2A D . The 
dash-dotted curve corresponds to the strongly collisional (SC) limit for a 
nonabsorbing particle. The dashed curve shows the potential for an absorbing 
particle in the same highly collisional limit. Both last curves are calculated 
using Eq. (14} for the same parameters as for the CL case but with T e = Tj, 
a/A D = 0.02, £i/X D = 0.1. 



by the presence of the particle. Equation (13\ is accurate to 
0(M T /r 5 ) + 0(M T /r 3 ) at M T -> 0, r -> oo. This result 
shows: 

• At sufficiently large distances the potential has the r~ 3 - 
dependence, which is in agreement with the inverse 
third power law of screening in anisotropic collisionless 
plasmas lETTl . 

• The resulting electric interaction between the particles 
is non-reciprocal {actio ^ reactio) since </>(r) ^ </)(— r). 
For instance, if two grains are aligned along the flow then 
the grain located downstream may experience attraction 
whereas the grain located upstream will be always re- 
pelled. 

• The potential in the direction perpendicular to the flow 
does not have an attractive part. This is contrast to 
the case of a shifted Maxwellian distribution (see, e.g., 
Eq. (19) of Ref. (62), Fig. 2 of Ref. ED, or Fig. 3 of 
Ref. [64]) and shows the importance of accounting for 
the non-Maxwellian form of the ion distribution. 

The wake effect has been proposed to be used to design 
binary interactions in complex plasmas by applying electric 
fields of various polarizations 11601 . The idea is to apply an 
electric field oscillating with a frequency which is (i) much 
faster than that characterizing the particle dynamics and (ii) 
much slower than the ion plasma frequency and the ion-neutral 
collision frequency. In this case the particles do not react to 
the field, whereas ions react instantaneously and the resulting 
interaction between the particles is determined by the time- 
averaged distribution of the electric potential in the screening 
clouds. Here the most interesting case is the case of spherical 
polarization where the vector of the electric field rotates in 
such a way that its directions are isotropically distributed in 



3D space but its absolute value remains constant. In this case 
the resulting potential is isotropic but the second and third 



terms in Eq. (13 1 are averaged out so that the effect of the 
field is in the rest term 0(M T /r 5 ) + 0(M T /r 3 ). The resulting 
potential was investigated numerically in Ref. [60 1 for finite 
Mt and found to have an attractive part. 

The effect of ion absorption on the particle surface has been 
so far investigated only for subthermal ion flows in highly 
collisional plasmas using the hydrodynamic approach 11651 . In 
this regime the absorption-induced ion rarefication behind the 
particle can overcome the effect of ion focusing and a negative 
space charge region develops downstream from the particle. 
This can have important consequences for particle motion in 
low-ionized plasmas ll66l and, therefore, let us discuss this 
issue in some more detail. 

A stationary negatively charged spherical point-like particle 
which is immersed in a quasineutral highly collisional plasma 
is considered. The electric field is sufficiently weak so that the 
ions are drifting with subthermal velocity while electrons form 
stationary background. Plasma absorption occurs on the grain 
surface, i.e. it acts as a plasma sink. Ionization/recombination 
processes in the vicinity of the grain are neglected. The 
potential distribution is obtained by solving Poisson equation 
coupled to the hydrodynamic equations for ions and Boltz- 
mann equation for electrons using the standard linearization 
technique under the assumption Mt <C kj^ii <C 1. Further, 
using the known asymptotic expression for the ion flux to 
an infinitesimally small grain (a <C Ad) in the continuum 
limit (li < a) E6), EH, lEl the expression for the potential 
downstream from the particle (6 = 0) can be written in the 
form 



exp(-rfc D ) 



(r) = Q 



2 - (r 2 k% + 2rk u + 2) cxp(-rfc D ) 



1 - 



2k 2 



k 2 



M T k 2 m 



fcpA 1 - exp(-rfcp) 
ko J r 

. 2 - [(! - klJk^)- 1 ^^ + 2rfc D + 2] exp(-rfc D 



(14) 

Let us briefly analyze the structure of the Eq. ( |T4) . The first 
term is the usual isotropic Debye-Huckel potential and the 
second term represents the anisotropic part of the potential 
behind a nonabsorbing particle. The third and fourth terms 
represent, respectively, the isotropic and anisotropic parts 
of the electric potential associated with the effect of ion 
absorption. In the absence of ion absorption the long-range 
asymptote of the electric potential can be written in the form 



<j)Ki ~{Quv/u pi r )(kr>i/kn) 



(15) 



Except of the factor (/cDj/^D) 4 this expression is identical to 
that obtained in Ref. l68l for the potential behind a slowly 
moving non-absorbing test charge in collisional plasma. The 
difference appears because in Ref. [68 1 only electron screening 
was taken into account. Thus, in contrast to the collisionless 
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regime, the potential of a non-absorbing particle in collisional 
plasmas exhibits oc r~ 2 long-range decay. In the presence of 
absorption the long-range asymptote of the anisotropic part of 
the potential can be written as 

(f> « -2(Q^/^ r 2 )(A :D4 /fc D ) 4 (fc De /A :D ) 2 . 

Thus, the effect of absorption induced ion rarefication com- 
petes with the effect of ion focusing and the amplitude of 
the anisotropic part of potential decreases, provided T e > Tj. 
More important is the contribution from the isotropic part as- 
sociated with absorption. It is always negative (for a negatively 
charged particle) and can be written as 

<t>*(Q/r)(km/kv) 2 , 

This contribution is most often dominant and completely 
determines the long-range behavior of the potential 1 651 . 

Figure [2] shows an example of calculating the electric 
potential downstream from the particle for the situations 
discussed above. The solid curve corresponds to the non- 
absorbing particle in collisionless plasma [Eq. ( 13 l]. Ions are 



focused downstream from the particle and a positive space 
charge region exists. The dash-dotted curve shows the potential 
distribution for a non-absorbing particle in highly collisional 



plasmas [first two terms of Eq. (14i]. Here collisions enhance 



ion focusing 11691 and the positive space charge region grows 
considerably. However, when absorption is included using Eq. 
( |14) , the absorption induced ion rarefication plays a dominant 
role and the potential is negative (dashed curve), in contrast 
to the previous cases. 

B. Suprathermal ion drifts 

In the highly suprathermal regime Mt 3> 1 the potential 
can be found analytically using the kinetic approach with the 
realistic constant-mean-free-path collision term |70|: 



x K a 



exv[it{r\\/ti)] 
'tx / t 2 + (VA D , cff ) 2 X(^ 



(16) 



Here rj_ is the distance in the plane perpendicular to the flow, 
rii is the distance along the flow (rn > and ni < along and 
against the flow, respectively), Ad, off = [eEti/ {^itne 2 )] 1 ! 2 is 
the effective Debye length, E is the field which drives the flow 
[it is related to the flow velocity via u = y/2eEii/(nmi)], 
Kq is the zero-order modified Bessel function of the second 
kind CD, 



x(t) = i - yT 

no - 2VT 



it 



it 



it, 
l 



da 



1 



(17) 



, [l + it(l -a 2 )] 2 it(l + it)' 

and the square roots must be taken with positive real part; the 
electrons are considered as a homogeneous background which 
is not perturbed by the presence of particle. Expression (16i 
works at all but very small angles with respect to the flow 
since it diverges logarithmically at r± — >• for rn > due 




Fig. 3. Contour plot of the potential around a small charged dust particle 
in a weakly ionized plasma with suprathermal ion drift driven by an electric 
field [Eq. J 1 61] for Xu.es/^i = 0.5. The charge is in the center and the drift 
direction is trom the left to the right. The distance and the potential are in 
units of Apj e ff and Q/Apj e ff> respectively. 



to the neglect of the thermal spread of neutral velocities. The 
contor plot of potential ( [T6| is shown in Fig. [3] 

The derived potential applies to the screening of charged 
particles in plasma presheaths where the ion flow velocity is 
in between the thermal velocity of neutrals and the Bohm 
velocity. This potential has been shown in Ref. COll to be in 
agreement with measurements of Konopka et al. l72l . 

The potential ( 16 1 at large distances is l70l 



Q^T>,eB cos( 



3/2 



(18) 



which is accurate to O 



„-3\ 



at r — > oo. This is an unscreened 



dipole-like field with the dipole moment |Q|A D e a/£i directed 
along the flow for Q < 0, although there is a difference 
from a pure dipole field due to the factor [2/(1 + cos 2 6) 



3/2 



The oc r -dependence in Eq. (18 1 is again different from 



the inverse third power law decay in collisionless anisotropic 
plasmas because of a finite collision length. Note that for 8 = 



Eq. (18 1 is identical to the result of hydrodynamic approach 
[Eq.J 15 1], provided we assume T e — > oo and u = eE/rriiV in 



Eq. ( 15 i. In both cases the dipole momentum is QeE / {miio* 



IV. Summary 

In this paper, we have summarized and discussed our recent 
results regarding the electric potential distribution around a 
small charged particle in weakly ionized plasmas. Different 
effects influencing the shape of the potential, including plasma 
absorption on the particle, ion-neutral collisions, plasma pro- 
duction and loss processes, and plasma anisotropy have been 
considered. The generic property of the results obtained so 
far is that sufficiently close to the particle the potential can 
be well approximated by the Debye-Hiickel (Yukawa) form. 
At longer distances, the potential usually exhibits power-law 
decay oc c/r n . The value of n (n = 1,2, or 3 in the 
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cases investigated), as well as the sign and magnitude of the 
parameter c, depend on the plasma and particle properties. 
These results have important consequences for plasma-particle 
and interparticle interactions and related phenomena, including 
phase transitions, phase diagrams, transport, waves, etc. 
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